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ABSTRACT: Starting from a careful analysis of the coupled Langevin equations for 
two interacting Brownian particles, we derive a method for extracting the binary friction 
tensor from the correlation function matrix of the instantaneous forces exerted by the bath 
particles on the fixed Brownian particles, and from the relaxation of the total momentum 
of the bath in a finite system. The general methodology, which circumvents the pitfalls 
associated with the inversion of the thermodynamic and long time limits, is applied to the 
case of two Brownian hard spheres in a bath of light spheres. 
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Ever since Green-Kubo (GK) formulae have been derived, expressing linear transport 
coefficients as time integrals of correlation functions of thermally fluctuating dynamical 
variables, it has been known that particular care must be exercised in evaluating these 
integrals from correlation functions for finite systems [1] . Strictly speaking, GK formulae 
yield non-zero results only provided the thermodynamic limit is taken before the upper 
limit in the GK integral is taken to infinity. To obtain sensible, non-zero results from 
the integration of correlation functions of finite systems, like those provided by Molecular 
Dynamics (MD) simulations of samples of N ~ 10 2 — 10 4 particles, a somewhat arbitrary 
upper cut-off r^y must be applied. For most transport coefficients, involving systems of 
identical or similar particles, the resulting values are not very sensitive to the precise value 
of 7~tv, since it is found that after a time of the order of the initial, fast relaxation of 
the system under study (typically a picosecond for dense fluids), the integral reaches a 
"plateau" value, which roughly coincides with the time beyond which the MD-generated 
correlation function drops below the noise level. 

However, the difficulty is less easily overcome when one considers the classic example 
of the friction coefficient ( exerted on a heavy Brownian particle by a bath of much lighter 
particles. ( is related to the time integral of the autocorrelation function (ACF) of the 
instantaneous force exerted by the bath particles on the Brownian particle. Recent MD 
simulations clearly show that no well-defined "plateau" value of the GK integrand is ob- 
served in systems invoving several hundred bath particles, so that the cut-off time becomes 
totally arbitrary [2,3]. In practice ( was determined from the relaxation of the total mo- 
mentum of the fluid, due to the collisions with a fixed Brownian particle of infinite mass 
M. In this letter, we will consider the case of two Brownian particles suspended in a bath 
of discrete light particles. We will show that finite size effects are even more subtle in this 
system and lead to spurious results for the computed friction tensor. We will then present 
a method to overcome these effects and obtain the correct friction tensor. This is a first 
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step in a statistical (first principles) approach to the hydrodynamic interactions in sus- 
pensions; such interactions are traditionnally derived from macroscopic hydrodynamics [4] 
which ignore the discrete nature of the bath, and are hence expected to fail at nanometric 
scales, like those explored by modern surface force machines [5]. 

The friction tensor ( relates the fluctuations of the forces acting on each suspended 
Brownian particle to the velocities of these particles : 



^a(t) = -i a V 1 (t)-gV 2 (t); a = 1,2 (i) 
where the 2x2 matrix £ of friction tensors is a function of the relative position R = 

=ab 

Rl — R2 of the two Brownian particles. The tensor occurs naturally in the two-particle 
Fokker-Planck equation, describing the dynamical evolution of two massive particles in a 
bath of much lighter particles [6,7,8]. The resulting microscopic expression for the friction 
tensor reads : 

1 f 00 

C =ab = k^fJ Q dT mKa] 0) ' 6nKb] r)) N|Ri,R 2 ) ( 2 ) 
where the notation (e<?|Ri, R2) refers to an equilibrium average over the fluid variables in 
the presence of two fixed Brownian particles located at Ri and R2; 5jF(R a ; t) = J^(R a ; t) — 
(•^ r (Ra))( eg |R 1 r 2 ) is the fluctuation of the force experienced by the Brownian particle at 
R a due to collisions with fluid particles. Due to the presence of the other Brownian 
particle, the average force experienced by particle a does not vanish, and is substracted 
from the instantaneous force, to yield the fluctuating force for a given configuration of the 
Brownian particles. Recently we have provided a rigorous derivation of the two-particle 
Fokker-Planck equation for a system of two Brownian spheres of diameter E in a fluid of 
N smaller spheres of diameter a [9]. Use of the multiple time-scale analysis, previously 
applied to derive the Fokker-Planck equation for a single Brownian particle [10,3], avoids 
any "ad hoc" assumptions concerning the separation of time scales. For hard spheres, the 
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fluctuating force reduces to the rate of transfer of momentum from the bath to Brownian 
particles in the course of instantaneous elastic collisions, as one might have intuitively 
expected [11]. The average force, on the other hand, which will henceforth be denoted by 
T a for brevity, may be identified with the familiar entropic depletion force acting between 
sterically stabilized colloidal particles [12]. 

The difficulties encountered in attempting to evaluate the friction tensor from eq. (2) 
for a finite system are immediately apparent if one notes that for a system of N fluid 
particles in the presence of two Brownian particles, the time derivative of the total fluid 
momentum may be identified with the sum of the forces acting on the two Brownian 
particles. It follows from this observation that for a finite system : 

i r°° 

£l fl + £ 2a = k^TJo dT ^ Ri;T ) +^( R 2^)] •<" r (R«;0)>(e 9 |R 1 ,R 9 ) 

^ t-00 (3) 

= k^fj MP-^(^;0)> (e(Z | Rl , R2 ) =o 

where P is the total fluid momentum. Since there is no reason why the diagonal and 
off-diagonal (or mutual) friction tensors should be exactly opposite, the spurious result (3) 
must be regarded as a consequence of the finite size of the system. 

In order to see how sensible results may be extracted from dynamical trajectories of 
finite systems, we consider the coupled evolution equations of the momenta P a of the two 
Brownian particles, in the form of generalized Langevin equations, namely : 



6=1,2 1/0 



(4) 



where is the matrix of memory functions, and SJ-"^(t) is the fluctuating "random" force, 
associated with the "fast" fluid variables; it satisfies the following constraints : 

(8F+(t)) = 0, (8F+(t)P b (t')) = 0, (8Ft{t)m')) = (5) 

Note that the mean forces T b depend on time through the slowly varying Brownian par- 
ticle positions. The clear separation of time scales associated with the heavy Brownian 
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particles (slow variables), and with the light fluid particles (fast variables), justifies the 
usual Markovian assumption, namely 

(6J*(t)6J*tf)) = Z£ atb 6(t-l!) = 2y™dr (^+(r)^+(0))} S(t-t') (6) 

which amounts to neglecting memory effects, so that the coupled generalized Langevin 
equations (4) reduce to the phenomenological local Langevin equations. Using matrix 
notations (e.g. V = (Pi, P2) etc.), they can be cast in the form 

p(t) = F(t) + M-V{t) + 5F + (t) (7) 

where the t-independent matrix }A_ is related to the ACF matrix of the random force via 
the fluctuation-dissipation theorem : 

& = g.(W)- 1 (8) 

Taking the Laplace transform of the Langevin eq. (7), and projecting onto 5jF(t = 0) leads 
to the following relation for the fluctuating force ACF : 

{6~F(s)6F) = |l+ = j • (5y + (s = 0)S^ + {0))=Ur 1 {s){ST + {s = 0)£F+(0)) (9) 

where the tilda denotes a Laplace transform and s is the variable conjugate to time. 

The inverse of the correlation matrix of the Brownian particle momenta may be cal- 
culated in the microcanonical ensemble (appropriate for MD simulations) for a system of 
2 Brownian particles of mass M and N bath particles of mass m, using the results of ref. 
[2]: 

where A = 2M+Nm • ^ ec l" (^) ^ s substituted into eq. (9) and the thermodynamic limit is 
taken before the Brownian limit M — > 00, the matrix A4 is found to vanish, and comparison 
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of relation (9) taken for s = 0, with the GK relation (2) yields the following expression for 
the friction tensor matrix : 

1 f°° 

i=k^Tj Q (^ + (0)6T + (r))dr (11) 

In other words, provided the limits are taken in the order specified above, the friction 
tensor is given equivalently in terms of the "bare" or "random" force ACF. For a finite 
system, the two estimates, say ( N and C^' + corresponding respectively to eqs. (2) and 
(11) differ. However, the spurious behaviour of the bare force ACF, embodied in the result 
(3) which is a direct consequence of the conservation of the total momentum of a finite 
system, is not expected to carry over to the random (or projected) force ACF, which are 
associated with the fast fluid variables. Consequently no singularity is expected when the 
order of limits (A — > oo, M — > oo) is inverted, and hence the finite and infinite system 
results for ( N > + should only differ by terms of order 1/N, i.e. : 

i N ^=i N ' + +°(Jr) ( i2 ) 

Discarding henceforth the superscript A for the tensor C + to simplify notations, we first 
note that symmetry considerations imply that the tensors ( + are diagonal and 

=ab 

£+ = = Ci (** + yy) + Cjf* (i3.a) 

^ 2 = ^ 1 = Ci(**+yy) + Cr (13 - b) 

where the indices s and m refer to "self" and "mutual" , and the Oz axis has been chosen 
along R = Ri — R2. 

Taking the Brownian limit M — > 00 for a system with finite A, it is easily found from 
eqs (6), (8) and (10) that the matrix U_~^(s = 0) appearing in eq (9) goes, when s — > 0, to 
the finite value 

{Ur 1 (s = 0)} ab = ^l-(28 ab -l) (14) 
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This non-vanishing result can be directly traced back to the conservation of total (fluid + 
Brownian) momentum. Consequently we derive from eq. (9) the following relation between 
the friction tensor matrices ( N and C + , valid for any finite N : 

C = k^f^ = °>^> = \ (£n - £ 2 V) • <■"* - v < 15 » 

In other words, keeping in mind eq. (12) : 



^■^■ te-&0 +o( ^i (16) 

=11 =12 2 \N ' 



1 



= ^ZT-^ (17) 



We thus have established a first relation between the true friction tensor and the tensor 
( N which may be calculated from the total force ACF of a finite system. From the 
small- s expansion of l£^(s) in eq. (9) one finds that the time integral of the force ACF 
(5jF(t)5jF(0)) relaxes exponentially towards its infinite time value. The relaxation times 
associated with the components of the fluctuating force perpendicular and parallel to the 
vector R are found in terms of the transverse and longitudinal components of the tensors 
(13), namely 

Nmk B T 

These relaxation times are seen to be proportional to the system size. However, one can 

easily check that the difference between the time-integral of the self and mutual force 

ACF's relaxes towards the difference of the corresponding friction matrices, — , on a 

much faster time scale, of the order of the correlation time of the random force. This fast 

relaxation is due to a compensation between the slow decays of the separate functions, and 

allows accurate estimates of the difference C + — C + from MD simulation, as illustrated in 

=11 =12 

Fig 1. 

At this stage, a second relation is needed to determine the self and mutual friction 
coefficients separately from MD simulations of a finite system. In a spirit similar to the case 
of a single Brownian particle [2,3], we now consider the relaxation of the total momentum of 
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N spheres in the presence of two infinitely massive (i.e. fixed) Brownian spheres. Summing 
the two Langevin equations (4) for Pi and P2 in the Markovian limit, noting that T\ = 
—Tli and averaging over a time interval intermediate between the short time scale of the 
fluctuating random force and the long time scale associated with the Brownian particles, 
we find for the evolution of the total momentum P = Pi + P2 : 

p(t) = _ 2 {c +C (18) 
v ' 1=11 =12 J mN v ' 

Due to the conservation of total momentum, (18) also holds for the momentum of the 
fluid P(£) = Nmv(t), where v(£) denotes the center-of-mass velocity relative to the fixed 
Brownian particles. We conclude that the components of the fluid momentum relax ex- 
ponentially, on time scales identical to those of the integrated ACF's, given in eq. (17). 
In practice, a MD calculation of the logarithm of the normalized fluid ACF : log F a (t) = 
should be a straight line with slope 2(^+0, + ^^/(Nrnkj^T), where 
a = {||,-L} is either a longitudinal or transverse component. An example of such a plot is 
shown in Fig. 2. The slopes thus provide the second relation between the longitudinal and 
transverse elements of the self and mutual friction tensors, which together with eq. (16), 
entirely determine the latter. 

The various coefficients depend on the distance R between the Brownian particles, and 
MD simulations must be carried out for different spacings of the fixed Brownian spheres. 
To illustrate the procedure, we quote the values of ( s ,a and Cm, a (a = {||>-L}) obtained 
for a size ratio of the Brownian and fluid spheres E/cr = 2, and a distance R/T, = 1. 
The packing fraction of the fluid is 77 = 0.246. The following values are extracted from 
plots as shown in Figs. 1 and 2 : C*y = C s ,\\/ mv mM = 12.1, C* ± = Cs,J mv mM = 6-7, 
= (m,\\/ mi/ mM = -0-3, Cm,± = Cm,±/ mv mM = -0.3, where v mM is the collision 
frequency between the bath and one Brownian particle. Contrary to the prediction of 
hydrodynamics, the value remains finite when the two Brownian particles are at contact. 

8 



(pJtyp a (o)) N 

i0 & Nmk B T 



In general our numerical results indicate strong deviations on these scales from the ill- 
dependence of the friction coefficient from hydrodynamics. More complete report, for 
various size ratios E/<r and distances R will be published elsewhere, and confronted with 
the predictions of macroscopic hydroynamics. 
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FIGURE CAPTIONS 



Figure 1 : Time-dependence of the reduced self (upper curve) and mutual (lower 
curve) friction coefficient, with appropriate sign (e a 5 = +1 if a = b and — 1 if a ^ b). 
Cfe(^) = Cab(^)/ mi, mM * s defined as the time-integral of the corresponding force ACF, up 
to time t (see eq. (2)), with v m M the fluid-Brownian particle collision frequency. The half 
difference between the two previous estimates (open circles) , is shown to converge towards 
the same asymptotic value on a much faster time-scale, of the order of a few collision time 
inside the fluid. 

Figure 2 : Logarithmic plot of the normalized ACF, F a (t), of the total fluid mo- 
mentum versus reduced time v m M^- The solid (resp. dashed) line corresponds to the 
component of the momemtum ACF perpendicular (resp. parallel) to R = R4 — R2. 
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